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Abstract 


The  model  M/G/l  is  modified  in  the  following  manner.  The  server  is  idled  when  he 
runs  out  of  customers  and  resumes  serving  when  the  (m+l)st  customers  arrives,  where 
m  £  0;  m+1  is  referred  to  as  the  ‘‘initial  quorum”' and  m  is  the  maximal  queue  size 
(and  system  size)  while  the  server  idles  or  vacations.  The  modified  model  is  designated 
as  M/G/l (m).  When  m=0  we  have  the  regular  M/G/l. 

In  Section  1  we  derive  the  omni-equations  for  the  backlog  process  B  (equation 
(1.17))  and  for  the  delay  process  w  (equation  (1.24))  by  exploiting  the  simple 
relation  between  B  and  w  in  M/G/l(m),  a  relation  which  results  in  equation  (1.12). 

In  Section  2  we  derive  several  composition  properties  for  the  backlog  and  for  the 
delay.  In  particular,  equation  (2.8)  says  that  the  backlog  is  distributed  like  the  sum  of 
the  backlog  in  M/G/l  and  of  a  random  variable  which  depends  on  the  service  time; 
and  equation  (2.9)  says  that  ihe  delay  is  distributed  like  the  sum  of  the  delay  in  M/G/l 
and  of  a  random  variable  which  depends  on  the  service  time  and  on  the  interarrival 
time.  -  .  .  ;  A  ~  ,,  "•  'r  .  .?  ■■  <  r  .  /'  >■ 

Notation 

A  =  interarrival  interval 

A.  =  arrival  rate  =  1/A 
x  =  service  duration 

M-  =  i/x 

p  =  \/p 

r  =  residual  time  (“residue”)  of  x;  x  and  r  are  related  through  the  omni-equation 
E<t>(x)  -  4>(0)  =  ExE4>'(r);  cf.  Krakowski  (Sept.  1984) 
w  =  delay  encountered  by  a  true  or  virtual  customer;  in  regular  M/G/l  we  have  w=B 
B  =  backlog  (unfinished  work)  due  to  incumbent  customers;  in  regular  M/G/l  we  have 
w  =  B 

w.  =  delay  when  server  works 

B,  =  backlog  when  server  works 

m  =  maximal  size  of  the  queue  when  the  server  idles;  m+1  is  the  “initial  quorum” 

N  =  size  of  the  system 
n  =  size  of  the  queue 
P,  =  Pr(N=j) 


Qj  =  Pr(n=j  and  server  idles)  ,  0  <  j  <  m  ;  Qq  =  PQ 
P,  =  Pr(server  works) 

X,  Y,  and  Z  are  defined  by  (1.5),  (1.6)  and  (1.23) 

4>  (  )  is  an  arbitrary  well-behaved  function  of  its  argument(s);  cf.  Krakowski  (1984, 
1985);  polynomials  and  their  limits  are  well-behaved  if  their  expectation  is  finite; 
so  is  the  step  function;  E4>(z),  the  expectation  of  4>(z),  is  called  the 
omni-transform  of  z. 

The  omni-convention  calls  for  taking  the  expectation  of  each-side  of  an 
omni-equation  without  explicitly  showing  the  expectation  operator;  thus  4>(y)  =  4>(z) 
stands  for  E4>(y)  =  E<j>(z) 

A  free  copy  of  a  random  variable  z  is  a  random  variable  distributed  like  z  but 
independent  of  z  and  of  any  other  variable  within  the  scope  of  the  same  expectation 
operator  (the  expectation  may  be  implied  by  the  omni-convention) 

j*z  =  zl  ++  Zj  where  j  is  a  positive  integer  and  the  zy  are  free  copies  of  z  (the 
“generic”  random  variable);  0*z  =  0 


Section  1.  The  Model  M/G/l  With  Initial  Quorum 


The  regular  model  M/G/l  is  modified  as  follows.  The  server  is  idled  (or  goes  off 
on  vacation)  each  time  he  runs  out  of  customers,  and  resumes  serving  at  the  instant  of 
the  arrival  of  the  (m+l)-th  customer,  where  m>0.  Thus  m+1  constitutes  an  initial 
quorum  and  m  is  the  maximal  queue  size  while  the  server  idles;  we  can  think  of  m 
as  a  “limbo.”  We  will  designate  this  model  as  M/G/l/(m).  Thus,  when  m  =  0  we  have 
the  regular  M/G/l. 

In  a  prior  report  (Krakowski,  1986)  we  dealt  with  the  process  “queue  size”  for 
M/G/l (m).  (Note:  In  that  prior  report  m  stands  for  the  current  m+1.  The  change 
simplifies  somewhat  the  typography  of  many  formulas.)  We  have  shown,  in  particular, 
that  for  each  0  <  j  <  m 

Qj  =  (l-p)/(m+l)  and  P.  =  p  =  X/jx  (1.1) 

where  Qj  =  Pr(n=j;  server  idles)  and  P.  =  Pr(server  works); 
and  that 


<t>  (n)  =  - — 2.  [4,  (o)  ++  4)  (m)]  +  p  4>  (n+h) 

m+l 


(1.2) 


where  h  =  number  of  arrivals  during  the  residue  of  x,  the  (generic)  service  time.  To 
make  this  report  readable  without  recourse  to  the  report  quoted  above  we  rederive 
(1.1).  (We  will  not  need  (1.2)  in  this  report.) 

Consider  the  transition  diagram  for  M/G/l/(3): 


server  works:  R  states 


server  idles:  Q  states 


The  flow  balance  for  the  state  “N=l;  server  idles”  is  XQ0  =  XQj  which  says  that  the 

frequency  of  entering  this  state  equals  the  frequency  of  exiting  it.  The  flow  balance  for 
the  state  “N=2;  server  idles”  is  XQj  =  XQ2;  and  for  the  state  “N=3;  server  idles”  the 

flow  balance  is  XQ2  =  XQ3.  Hence  we  have  Q0=Qi=Q2=Q3-  Clearly,  for 
arbitrary  m  we  have  Q0  =  Q-;  0  <  j  <  m.  The  value  of  P.  follows  from  the  global 
balance  of  arrivals  and  departures:  X  =  jj.P..  Thus  we  have  demonstrated  (1.1). 


In  the  model  M/G/l/(m)  the  treatment  of  the  delay,  the  theme  of  the  current  report, 
is  more  complex  than  the  treatment  of  the  queue-size.  It  is  also  more  complex  and 
more  subtle  than  the  analysis  of  the  delay  for  the  regular  M/G/l.  The  instinctive 
approach  in  analyzing  the  delay,  in  M/G/l/(m)  as  well  as  in  M/G/l,  is  to  establish  a 
balance  equation  for  the  virtual  delay.  This  has  been  a  simple  task  in  M/G/l,  mainly, 
in  our  opinion,  because  in  this  model  the  process  “virtual  delay”  and  the  process 
“virtual  backlog”  are  identical.  In  M/G/l/(m)  when  m>l  the  two  processes  differ.  It 
turns  out  that  the  analysis  of  M/G/l/(m)  is  much  simplified,  and  more  insightful,  if 
one  analyzes  the  backlog  process  B  along  with  the  delay  process  w;  the  backlog  B 
is  the  aggregate  remaining  service  time  for  all  incumbent  customers.  (Still,  the  direct 
“instinctive”  approach  of  balancing  the  virtual  delay,  or  rather,  in  the  omni-fashion,  an 
arbitrary  function  of  this  delay,  is  a  viable  approach.)  And,  of  course,  the  process  B 
may  be  of  interest  in  its  own  right  and  not  as  a  mere  stepping  stone  to  analyzing  w. 
Note  The  actual  delay  experienced  by  a  true  customer  has  the  same  distribution  as 
the  virtual  delay;  cf.  Wolff  (1982).  This  follows  from  the  fact  that  a  true 
customer  and  a  random  observer  see,  stochastically  speaking,  the  same 
picture  when  the  arrival  source  is  poissonian.  Using  the  same  symbol  w  for 
true  and  virtual  delays  should  cause  no  ambiguity  in  the  context  of  this 
report. 

The  Static  Balance  A  “random”  observer  will  find,  with  probability  Q0,  that  B  =  0  and 
N  =  0;  he  will  find,  with  probability  Qlt  that  N  =  1  and  the  server  idles  and  B  =  x; 
and  more  generally  he  will  find,  with  probability  Qj  where  0  <  j  <  m,  that  N  =  j  and 
the  server  idles  B  =  j*x;  and  he  will  find  with  probability  P,  that  B  =  B,.  Thus,  since 
each  Qj  =  (l-p)/(m+l)  and  P,  =  p,  we  have 

4>  (B)  =  (4>(0*x)  ++4>(m*x)]  +  P4>  (B.) 

m+ 1 

When,  e.g.,  m  =  3  equation  (1.3)  becomes 

4>(B)=  — 2-  [4>(0)  +  4>(x)  +  4)(2*x)-t-4>(3*x)]  +  p4>(B.)  (1.3a) 

4 

To  find  the  omni-equations  for  the  delay  w  which  correspond  to  (1.3)  and  (1.3a)  we 
note  that  a  virtual  customer  who  finds,  with  probability  Qj,  j  (true)  customers  in  the 

limbo  must  await  m  -  j  additional  customers  to  reach  the  quorum  of  m  +  1  (this 
quorum  includes  the  virtual  customer)  and  then  must  wait  until  the  j  senior  customers 
are  served.  The  entire  delay  of  our  virtual  customer  is  therefore  (m-j)*A  +  j*x,  where 


4 


0  <  j  <  m;  thus  the  delay  of  a  cusomter  who  finds  the  server  idle  is  composed  of  m 
intervals.  (When  m  =  0  there  is  no  delay  when  the  server  idles,  this  being  the  regular 
M/G/l.)  Hence, 

<J> (w)  =  [(j>(m*A)  +  4>((m-l)*A+x)  ++4>(m*x)]  +  P<j>  (w, )  (1-4) 

m+ 1 

When  m  =  3  equation  (1.4)  becomes 

4>  (w)  =  [4>(3 *A>  +  4>(2*A+x)  +  <{)(A+2*x)  +  4>(3*x)]  (l-4a 


+  p4>(w.) 

where  w,  is  the  delay  conditioned  upon  the  server  working. 
We  simplify  the  typography  of  (1.4)  and  (1.4)  by  defining 

■MX)  a-!-  g  <K j •  x) 
mtl  ,fo 

and 


4>(Y)  i-L- 

m+l 

Then  (1.3)  and  (1.4) 


§  4>((m-j)*A+j  *x) 
i=o 

become 


(1.5) 


(1.6) 


4>  (B)  =  ( 1-p)  4>(X)  +  P4*  (B, )  (1.7) 

and 

4>(w)  =  (1-p)  <t>(Y)  +  p<j>(w.)  (1.8) 

A  key  observation  in  our  approach  is  that  in  the  model  M/G/l/(m)  the  virtual  delay 
and  the  backlog  observed  at  an  instant  when  the  server  works  are  one  and  the  same 
process;  hence 


B.  =  w, 

from  which  follows  the  weaker  statement  that 


(1.9) 


4>(B.)  =  4>(w. )  (1.10) 

which  merely  says  that  B,  and  w,  are  identically  distributed.  Hence,  from  (1.8  and 
(1.10)  it  follows  that 

4>(w)  =  (1-p)  4>(Y)  +  p4>(B, )  d  ll) 
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It  follows  from  (1.7)  and  (1.11)  that 


<J>(w)  =  4>(B)  +  (1-p)  (4>(Y)  -4>(X)]  (1.12) 

Thus  if  we  have  an  omni-equation  for  B  we  can  transpose  it  with  the  aid  of  (1.12) 
into  an  omni-equation  for  w.  We  now  derive  such  an  equation  as  balancing  condition 
for  an  arbitrary  function  of  B. 

The  backlog  B  jumps  to  B  +  x  with  each  arrival,  these  arrivals  being  of 
frequency  k.  (Departures  of  serviced  customers  do  not  affect  the  balance  since  in  our 
model  we  assume  that  the  value  of  a  customer’s  service  time  is  revealed  at  the  instant 
of  his  arrival.)  On  the  other  hand,  while  the  server  works,  with  probability  P.,  he 

keeps  on  reducing  the  backlog  now  denoted  by  B,  at  the  rate  dB.  =  dt.  Hence  the 
omni-balance  for  B  is 

M<KB  +  X)-<KB)]  =  P.<j>'(B.)  (1.13) 

According  to  the  shifted  renewal  equation  (Krakowski,  September  1984) 

<t>(B  +  x)  -  4>(B)  =  x<j>'(B+r)  (1.14) 

where  r  is  the  residue  of  x;  the  differentiation  is  with  respect  to  the  entire 
argument.  From  (1.13)  and  (1.14)  we  have 

\x4>'(B+r)  =  P.4>'(B.) 
which  integrated  typographically  is 

p<J)(B+r)  =  P,<J)(B.)  (1.15) 

Putting  4>(  )  =  1  w'e  rederive  P,  =  p  and  (1.15)  becomes 

4>(B.)  =  4>(B+r)  (1.16) 

From  (1.17)  and  (1.16)  we  have 

4>(B)  =  ( 1-p)  4>(X)  +  P4>(B+r )  (1.17) 

the  sought  for  omni-equation  for  the  backlog  B.  From  (1.17)  we  can  get,  recursively, 
successive  moments  of  B  and  a  convolution  equation  for  the  distribution  of  B.  And 
in  turn,  with  the  aid  of  (1.12),  we  can  get  the  successive  moments  of  w  and  a 
convolution  equation  for  the  distribution  of  w.  Of  course,  all  needed  moments  and 
distributions  of  X  and  Y  and  Z  have  to  be  derived  as  a  side  exercise. 
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Example:  Find  B  and  w 
From  (1.15)  we  have 


X  =  — —  [x  +  +  mx  ]  =  JSiSliliS  =  1  mx  (1.18) 

rn+l  1  1  2  (m+1)  2  v  ' 

and  from  (1.6)  we  have,  similarly, 

Y=jmA  +  jmx  =  >m(Ui)  (1.19) 

From  (1.17)  and  (1.18)  we  have 


B='mx  +  _PL  (1.20) 

‘  i-P 

From  (1.12),  (1.18)  and  (1.19)  we  have 

w  =  B  +  (l-p)(Y- X)  =  +  (l-p)mA 

=  fl  +  rf  C-2') 

Thus,  w  is  composed  of  m/A.  and  the  expected  waiting  time  for  a  regular  M/G/l. 
This  is  not  incidental  and  we  will  take  up  the  problem  of  composition  properties  in  the 
model  M/G/l (m)  in  Section  2. 

To  complete  the  current  section  we  derive  the  omni-equation  for  w.  From  (1.12) 
and  (1.17)  we  have 

4>(w)  =  (l-p) (4>(Y) -p[4)(Y+r)  - <t>(X+r)]}  +p4>(w+r)  (1.22) 

A  comparison  of  equation  (1.22)  with  the  quite  simpler  (1.17)— even  if  we  replace 
4>(X)  by  its  explicit  representation  in  (1.5)— will  give  the  reader  an  idea  about  the 
more  extensive  need  for  algebraic  maneuvering  in  a  direct  derivation  of  (1.22)  from 
the  omni-balance  of  the  virtual  delay. 

The  formal  appearance  of  (1.21)  is  simplified  by  defining 


4>  (Z)  d  4>(Y)  -  p  [4>(Y+r)  -  <j)(X+r)] 
From  (1.22)  and  (1.23)  we  have 

4>(w)  =  (1- p)  <(>(Z)  +  p<J>(w+r) 


(1-23) 


(1-24) 


It  is  easy  to  verify  that  ( 1 .24)  leads  to  w  as  given  in  (1.21). 


Section  2:  Composition  Theorems 

When  m  =  0,  i.e.  when  we  consider  the  regular  M/G/l/(0)  alias  M/G/l,  we  have 
according  to  (1.5),  (1.6)  and  (1.23) 


4>(Xo)=0,  <MYo)=0,  4>(Z0)  =  0,  (2.1) 

We  will  use  in  this  section,  as  we  have  just  done  in  (2.1),  a  subscript  serving  as  a 
reminder  of  the  size  of  m  in  the  model.  This  is  not  a  necessity,  just  a  convenience; 
most  equations  in  this  section  relate  a  general  M/G/l/(m)  model  to  the  regular 
M/G/l/(0)  so  that  the  reminder-subscripts  will  identify  the  models.  Thus  we  will  write 
Xm,  Ym,  Zm,  Bm,  and  wm,  in  place  of  X,  Y,  Z,  B  and  w.  Equations 
(1.17)  and  (1.22)  for  m  =  0  are  written  now,  taking  into  account  (2.1), 

4>(B0)  =  (1-pH(0)  +  p<l)(B0+r)  (2.2) 

and 

4>(w0)  =  (1- p)4>(0)  +  p4>(w0+r)  (2.3) 

Equations  (2.2)  and  (2.3)  imply,  as  can  be  shown,  that 

4>  (B0)  =  4>  (w0 )  (2.4) 

which  we  know  from  the  structure  of  the  regular  M/G/l. 

Let  us  now  shift  (2.2)  by  Xm  thus  obtaining 

4>(B0  +  Xm)  =  (l-p)4>(Xm)  +  pcj)(B0+r  +  Xm)  (2.5) 


From  (2.5)  and  (1.17)  we  haave 

4>(Bm)-p4>(Bm+r)  =  <j)(Bo+Xm)-  p4>(Bo  +  r  +  Xm)  (2.6) 

Defining 

^(Bm)  |  4>(Bm)-p4>(Bm  +  r)  (2.7) 

we  write  (2.6)  as 

■KBm)  =  >V(B„  +  Xm)  1 


or,  reverting  to  <J>(  )  in  place  of  ),  as 

4>(Bm)  =  4>(Bo  +  xm  )  (2.8) 

Equation  (2.8)  states  that  the  backlog  Bm  in  M/G/l/(m)  is  distributed  like  the 
convolution  of  B0  ,  the  backlog  in  M/G/l,  and  Xm. 

From  (2.3)  and  (1.24)  we  derive  in  like  manner 

4>(wm)  =  4>(w0  +  Zm  )  (2.9) 

thus  showing  that  the  delay  in  M/G/l/(m)  is  distributed  like  the  convolution  of  w0  , 
the  delay  in  regular  M/G/l,  and  Zm 

The  composition  equations  (2.8)  and  (2.9)  should  be  especially  useful 
computationally  when  the  distribution  or  moments  of  the  delay  w0  are  known  and 
when  m  is  treated  as  a  parameter.  Of  course,  the  distribution,  or  moments,  of  Zm 
have  to  be  evaluated  as  a  side  exercise.  Similar  remarks  hold  for  the  distribution  and 
moments  of  B  should  the  backlog  be  a  process  of  interest.  The  labor  involved  may 
be  formidable  but  in  some  practical  applications  it  might  be  justified  economically. 

Another  composition  property  for  the  delay  wm  is  obtained  as  follows.  Providing 
(1.12)  with  m-subscripts  we  have 

«M'vm)  =  ( l  -P)  [4>  (Ym )  -  4>(Xm )]  +  4>(Bm ) 

Since  in  the  regular  M/G/l,  i.e.  in  M/G/l /(0),  B0 
form 

<MBm)  =  4>(wo  +  Xm)  (2.11) 

From  (2.10)  and  (2.11)  we  have 

4>(wm)  =  (l-p)I4>(Ym)-4>(Xm)l+<t>(w0+Xm)  (2.12) 

Equation  (2.12)  appears  simpler  numerically  than  (2.9),  especially  for  the  computation 
of  moments. 


(2-10) 

=  w0  we  can  write  (2.8)  in  the 
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